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ABSTRACT 

We use the direct Fourier method to calculate t he redshift-space power spectrum of the 
maxBCG cluster catalog (iKoester et al.l [20071 ) - currently by far the largest existing 
galaxy cluster sample. The total number of clusters used in our analysis is 12, 616. After 
accounting for the radial smearing effect caused by photometric redshift errors and 
also introducing a simple treatment for the nonlinear effects, we show that currently 
favored low matter density "concordance" ACDM cosmology provides a very good fit 
to the estimated power. Thanks to the large volume (~ OAh~^ Gpc'^), high clustering 
amplitude (linear effective bias parameter beS ~ 3 x (O.SS/crs)), and sufficiently high 
sampling density ('^ 3 x 10^^ ft."^ Mpc^"^) the recovered power spectrum has high enough 
signal to noise to allow us to find weak evidence 2a CL) for the acoustic features. 
These results are encouraging in light of the several propose d large cluster survey s. In 
case we use the photometric redshift errors as suggested in iKoester et al.l (|2007[ ) we 
are left with the excess large-scale power which has previously been noticed by several 
other authors. 

Key v^rords: cosmology: observations - large-scale structure of Universe - galaxies: 
clusters: general - methods: statistical 



1 INTRODUCTION 

Since the flight of the Cobe satellite in the early 1990's 
the field of observational cosmology has witnessed ex- 
tremely rapid progress which has culminated in the estab- 
lishment of the St andard Model for cosmology - the "con - 
cordance" model (Bahcall ct al."l999l: [s"pergel et al.ll2003h . 
This progress has been largely driven by the precise mea- 
surements of the angular temperature fluctuations of the 
Cosmic Microwave Backgro und (CMB) llHa nanv et al.ll2000l: 
Ncttcrficld ct al. 2002; Be nnett et all [200 3: Grain ge et al] 
20031 : iPearson et al.l 120031 ). However, in order to break sev- 
eral degeneracies inherent in the CMB measurements one 
has to complement this data with other sources of informa- 
tion, such as the mea surement s of the SNe la luminosity 
distances (jRiess et al. [[1998; Pc rlmutter et al. I [1999) or with 
the measurements of the large-scale structure (LSS) of the 
Universe as traced by galaxies or galaxy clusters. 

The simplest descriptor one can extract from the LSS 
measurements is the matter power spectrum. The broad- 
band shape of this spectrum is sensitive to the shape param- 
eter r = Qmh, and thus is useful in helping to establish the 
currently favored low matter density "concordance" model. 



as well as the amount of baryons in relation to the total mat- 
ter fb — Qb/^m- Currently the two largest galaxy redshift 
surveys are the 2dF Galaxy Redshift Surve£[ (2dFGRS) 
and the Sloan Digital Sky SurvejQ (SDSS) with its latest 
data release 5 (DR5), providing redshifts to ~ 220,000 and 
~ 675, 000 galaxies, respectively. The SDSS galaxy sample 
consists of two broad classes: (i) the MAIN galaxy sample, 
reaching redshifts of 2: ~ 0.25; (ii) the Luminous Red Galaxy 
(LRG) sample covering redshift range z ~ 0.15 — 0.5. Some 
other characteristics of these surveys are listed in Table [T] 

In addition to the precise measurement of the broad- 
band shape of the power spectrum, SDSS LRG and 2dFGRS 
samples have proven useful in allowing the detection of theo- 
retically predicted oscillatory features in the spectrum. The 
source of these spectral fluctuations is the same as that giv- 
ing rise to the prominent peak structure in the CMB angular 
power spectrum - acoustic oscillations in the tightly cou- 
pled baryon-photon fluid pr ior to the epoch o f recombina- 
tion ( Sunyaev fc Zeldovich 1970l : [Peebles fc Yul[T970 ). How- 
ever, the relative level of fluctuations in the matter power 
spectrum is strongly reduced compared to the correspond- 
ing fluctuations in the CMB angular spectrum, owing to the 
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Table 1. Some characteristics of the 2dFGRS, SDSS DR5, and maxBCG samples. 



Survey 


Number of 


Sky area 


Redshift 


Comoving volume 


Effective bias parameter wrt to 




objects 


{deg2) 


coverage 


(h-^Gpc^*) 


the model with erg = 0.85 


2dFGRS 


~ 220, 000 


~ 1200 


£ 0.3 


~ 0.07 


~ 1.0 


SDSS DR5 


~ 675, 000 


~ 5700 


MAIN: < 0.25 


MAIN: ~ 0.2 


MAIN: ~ 1.2 








LRG: ~'^.15 - 0.5 


LRG: ~ 1.3 


LRG: ~ 2.0 



maxBCG ~ 13, 800 ~ 7000 0.1 - 0.3 ~ 0.4 ~ 3 



fact that most of the matter in the Universe is provided by 
the cold dark matter (CDM) that does not participate in 
acoustic oscillations. After recombination, or more precisely 
after the end of the so-called "drag-epoch" which happens 
at redshifts of a few hundred, as the baryons are released 
from the supporting pressure of the photon gas they start 
to fall back to the CDM potential wells that have started to 
grow already since the matter-radiation equality. Since the 
amount of baryons in the total matter budget is not com- 
pletely negligible - baryons make up roughly one fifth - the 
acoustic structure imprinted onto the baryonic component 
also gets transferred to the spatial distribution of the CDM. 
After baryonic and CDM components have relaxed one ends 
up with ~ 5% relative fluctuations in the total mater power 
spectrum. 

In order to observe these relatively small fluctuations 
one needs to have: 

(i) Large redshift surveys, to reduce cosmic variance; 

(ii) High enough spatial sampling density, to decrease dis- 
creteness noise. 

Both of these criteria are currently best met by the SDSS 
LRG sample, and indeed, at the beginning of 2005 the detec- 
tion of the acoustic peak in the spatial two-point co rrelation 
function was announced by lEisenstein et al.l |200^. At the 
same time the final power spectrum measurements of the 
2dFGRS also rev ealed the existence of the acoustic features 
l|Cole et al.ll2005l ). Baryonic acoustic oscillations (BAO) in 
the powe r spect r um of the SDSS LRG (DR4) sample were 
found by iHiitsil (l2006al ) and the corresponding cosmologi- 
cal parameter estimation was carried out in [Hiitsi ( 2006c|). 
Also, several other d e tectio ns of BAO have b e en rep orted: 
iPadmanabhan et al.l l|2006l ) an d iBlake et all (l2007f) used 
SDSS photometric LRG sample, iTegmark et al (|2006|) ana- 
lyzed spectroscopic LRG sample, and finally 'Percival et al.l 
1 20071 ) carried out a combined analysis of the SDSS MAIN 
and LRG spectroscopic galaxy samples. 

The usefulness of BAO arises from the fact that it can 
provide us with a standard ruler in the form of the size of 
the sound horizon at decoupling, enabling us to carry out a 
purely geometric cosmological test: comparing the apparent 
size of the ruler along and perpendicular to the line of sight 
with the physical size of the sound horizon (which can be 
well calibrated from the CMB data) one is able to find Hub- 
ble parameter H{z) and angular dia meter distance dA(z ) 
corresponding to the redshift z (e.g. IHu fc HaimanI |2003| ) . 
Obviously, having determined H{z) and dA{z) one can put 



constraints on the dark energy (DE) equation of state pa- 
rameter wde{z). In that respect radial modes, which enable 
us to find H{z), are more useful since H{z) is related to 
wde{z) through a single integration whereas dA{z), which 
is given by the angular modes, involves double integration 
over WDEiz). Thus one would certainly benefit from accu- 
rate redshift information. 

If one wishes to use the BAO signal as a precise standard 
ruler there are a few complications one has to overcome: 

(i) As the density field goes nonlinear couplings between 
Fourier modes will modify the BAO signal, leading to the 
damping of the oscillations; 

(ii) In order to fully exploit the information in 
galaxy /cluster surveys one has to understand the relation 
of these objects to the underlying matter density field, i.e. 
one has to have a realistic model for biasing. 

These two complications can be fully addressed only 
through costly N-body simulations. There are several pa- 
pers that have investigated the detectability and possible 
systematics of BAO extrac t ion u si ng the numeri c al sim - 
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significantly more work towards these directions is probably 
needed before we can harness the full power provided by the 
potentially very clean geometric tool in the form of BAO in 
the matter power spectrum. 

It is worth pointing out that actually one of the first 
possible hints for the existence of BAO in the matter power 
spectrum came from the a nalysis of the spat ial clustering 
of the AbeU/ACO cluster s (iMiller et al.ll200lM ^ The power 



spectrum measurement in lMiller et al.l ( 200ll ) used 637 clus- 
ters with richness R ^ 1. In this paper w e analyze the 
maxBCG cluster sample (jKoester et al.|[2007^ which is cur- 
rently by far the largest cluster catalog available, spanning a 
redshift range 2 = 0.1 — 0.3, containing 13, 823 objects, and 
covermg ~ 0.4/i~^Gpc^ of comoving volume (see Table [1] 
for comparison with 2dFGRS and SDSS DR5). The biggest 
advantage of using galaxy clusters instead of the typical, i.e. 
L»-galaxies, is the fact that their spatial clustering signal 
is strongly amplified with respect to the clustering of the 
underlying matter, and thus one can achieve the same sig- 
nal to noise measurement of the power spectrum by using 
correspondingly less number of objects. In Fig. [T] we have 
compared the effective volume of the maxBCG survey with 
the two other surveys that have lead to the detection of 



^ To be more precise again, at the end of the "drag-epoch" . In 
the "concordance" ACDM model the sound horizon at the "drag- 
epoch" is ~ 5% larger compared to the one at recombination. 



* Interesting features in the spatial two-point correlation func- 
tion of the AbeU/A CO cluster sample were also revealed by 
lEinasto et"all lll997l '). 
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Figure 1. The effective volume (see Eq. (|2]l) of tlie maxBCG 
cluster sample in comparison to the 2dFGRS and SDSS LRG 
DR,4 samples. 



BAO: 2dFGRS and SDSS LRG. Under the assumption of 
Gaussianity the power spectrum measurement errors AP 
are simply found as: AP/P = 2/Vcs/Vk, where Vk is the 
volume of the shell in fc-space over which the angular av- 
erage is taken, i.e. 14 = 47rfc^ Afc/(27r)"^. It is no surprise 
that the SDSS LRG sample is currently unbeatable. How- 
ever, we notice that maxBCG sample should provide better 
measurement of the power spectrum on scales larger than 
k ~ 0.1 /iMpc"^ compared to the 2dFGRS. The rapid fall of 
the effective volume of the maxBCG sample on smaller scales 
is caused by the photometr i c reds hift errors of ~ 0.01 as 
estimated bv lKoester et al.l (|2007h . 

Thus one would expect the maxBCG sample to reveal 
acoustic features. As our further analysis shows this indeed 
turns out to be the case, albeit with a relatively modest 
confidence level. 

Our paper is organized as follows. In Section 2 we give 
a brief description of the maxBCG catalog with the corre- 
sponding selection effects. The core part of this paper, which 
is devoted to the power spectrum analysis, is given in Section 
3. Finally, Section 4 contains our conclusions. 



DATA AND SURVEY SELECTION 
FUNCTION 



A maxBCG cluster catalog l|Koester et al.|[2007h is com- 
piled via the "r ed seque nce" clust er detec tion method 
llOstrander et al.| [l998: Gla dders fc Ye c 2000; BahcaU et all 
|20G3|) applied to the SDSS photometric data. The catalog 
contains 13, 823 galaxy clusters with velocity dispersions 
greater than ~ 400 km/s and covers ~ 7000 square de- 
grees of sky between redshifts 0.1 and 0.3. The photometric 
redshifts are estimated using the tight relation between the 
ridgeline color and redshift and are claimed to have accuracy 
Sz = \J {(zphoto — Zspcc)^) ~ 0.01 essentially independent of 



^ http:/ /umsdss. physics. lsa.umich.edu/catalogs/ 
maxbcg_public_catalog.dat 
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Figure 2. Angular distribution of 12,616 maxBCG clusters in 
SDSS survey coordinates (r), A). The 5x5 chessboard pattern of 
black/gray points shows the angular part of the division used for 
"jackknife" error analysis. The union of solid rectangles represents 
our reconstruction of the angular mask. 



redshift. For more details about the catalog we refer to the 
original source iKoester et al.l (|2007l ). 

In our power spectrum analysis we are going to neglect 
the three southern SDSS stripes, leaving us with 12, 616 
galaxy clusters over ~ 6800 square degrees of sky. With 
the redshift range Q.l < z < 0.3 this corresponds to 
~ 0.4 Gpc"^ of comoving cosmic volume. Neglecting the 
three narrow southern stripes helps us to achieve a better 
behaved survey window function with reduced sidelobes, and 
thus with reduced leakage of power. 

The angular distribution of the remaining clusters is 
given in Fig. (2] Here the angular coordinat es are plotted us- 
ing th e SDSS survey coordinates (??, A) fe.g. lStoughton et al.l 
l2002l l. The chessboard pattern of black/gray points in Fig. 
[2] represents the division of the sky into 5x5 angular re- 
gions used for "jackknife" error analysis. We also use three 
divisions along the redshift direction, making the total of 
3 X 5 X 5 = 75 regions, each containing approximately 168 
clusters. The union of solid rectangles in Fig. [5] represents 
our reconstruction of the survey angular mask. All the rect- 
angles are aligned with the SDSS imaging scan stripes. As 
the number density of galaxy clusters in the maxBCG sam- 
ple is rather high, one can determine relatively accurately 
the beginning, ending, and also possible gaps in the scan 
stripes. We call the angular mask obtained this way a "min- 
imal" mask in contrast to the "maximal" mask, which is 
built in the same manner except that each of the rectangles 
is extended by an amount corresponding to the mean cluster 
separation. We carry out our power spectrum analysis using 
both of these angular masks. 

In Fig. [3] we show the redshift distribution of maxBCG 
clusters. Here the solid smooth line corresponds to the cubic 
spline fit and the dashed line represents the expected num- 
ber of clusters per redshift bin = 0.2/30 for the volume- 
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Figure 3. Redshift histogram of the maxBCG sample. The sohd 
smooth hne corresponds to the cubic sphne fit and the dashed 
line represents the expected number of clusters per redshift bin 
for the volume-limited sample. 



limited sampl^fl assuming a spatially flat ACDM cosmology 
with flrn = 0.27. To convert this redshift distribution to 
the radial selection function n{r), i.e. the comoving number 
density of clusters at distance r, we fix the background cos- 
mology the low matter density model mentioned above. Our 
power spectrum calculations are done for both of the radial 
selection models shown in Fig. [S] 

For the full survey selection function n(r) we assume, 
as usual, that it can be factorized to the product of angular 
and radial parts, i.e. n(r) = n(f)n(r) (r = |r|, f = r/r). 
Here n{r) is the radial selection function and the angular 
selection part n(f) is assumed to take value 1 if the unit 
direction vector f is inside the survey mask and otherwise, 
i.e. we assume a simple binary angular mask. 



3 POWER SPECTRUM ANALYSIS 

We calculate the redshift-space power spectrum of the 
maxBCG catalo g using t he direct Fourier method as de- 
scribed in .Feldman et al. 1 ((1994 ) (FKP). Strictly speaking, 
power spectra determined in this way are so-called pseu- 
dospectra, meaning that the estimates derived are convolved 
with a survey window. Since in the case of the analyzed 
maxBCG sample the volume covered is very large, reach- 
ing OAh~^ Gpc'^, and also the survey volume has relatively 
large dimensions along all perpendicular directions, the cor- 
relations in the Fourier space caused by the survey window 
are rather compact. On intermediate scales and in the case 
the power spectrum binning is chosen wide enough, FKP es- 
timator gives a good approximation to the true underlying 
power. 

Instead of dir ect summation as presented in 
iFeldman et al.1 (|1994| ) we speed up the calculations using 



Fast Fourier Transforms (FFTs) . This gives rise to some ex- 
tra complications. As our density fiel d, which is built using 
the Triangular Shaped Cloud (TSC) (|Hocknev fc Eastwood 
,1988. ) mass assignment scheme, is now "restricted to 
live" on a regular grid with a finite cell size, we have to 
correct for the smoothing effect this has caused. Also, if 
our underlying density field contains spatial modes with 
higher frequency than our grid's Nyquist frequency, feNy, 
then these will be "folded back" into the frequency interval 
the grid can support, increasing power close to fcNy - 
the so-called aliasing effect. Quite often in the literature 
this aliasing effect is not properly accounted for, which 
is only fine in cases where the grid's Nyquist frequency 
used in the analysis is significantly higher than the spatial 
frequencies where one wants to measure the spectrum. We 
correct for the alias i ng eff ect using the iterative scheme 
as described in |jing| (j200fj) with a slight modification: we 
do not approximate the small-scale spectrum by a simple 
power law, but also allow for the possible running of the 
spectral index, i.e. the parametric shape of the power 
spectrum is taken to be a parabola in log-log. 

For the full details of our power spectrum calculation 
method along with several tests we refer the reader to iHiitsil 
l|2006bl lal). As a very brief summary, our spectrum determi- 
nation consists of the following steps: 

• Determination of the survey selection function, i.e. 
mean expected number density without any clustering n(r) 
(see Section [2}. This smooth field, with respect to which 
the fluctuations are measured, is modeled using a random 
catalog that contains 100 times more objects than the real 
maxBCG catalog, i.e. 1, 261, 600 in total; 

• Calculation of the overdensity field on a 512'^-grid using 
TSC mass assignment scheme; 

• FFT of the gridded overdensity field; 

• Calculation of the raw 3D power spectrum by taking 
the modulus squared from the output of the previous step; 

• Subtraction of the shot noise component from the raw 
3D spectrum; 

• Recovery of the angle averaged s pectrum usin g a mod- 
ified version of the iterative scheme of Ijingj l|2005h . 

The results of our power spectrum calculation along 
with fitted models are shown in Figs.|4]and[5] Here the error- 
bars AP are calculated as described in FKP, which assumes 
that the density fiel d follows Gaussian statistics (see also 
iTegmark erahU 19981 ): 



AP _ 

where the effective volume T4ff (see Fig. [l| is given by: 
n(r)P{k) 



l + n{r)P(k) 



dr , 



(1) 



(2) 



and Vk is the volume of the fc-space shell over which the 
angular average of the 3D spectrum is taken, i.e. 



Vk 



(27r)3 



(3) 



^ In iKoester et al.l ([SooJ) it is suggested that maxBCG sample 
is rather close to being volume-limited. 



In Eq. ([2]) Tmin and Tmax are comoving distances correspond- 
ing to the redshift bounds of the survey (0.1 and 0.3 in case 
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Figure 4. Upper panel: Filled circles with error-bars show the 
redshift-space power spectrum of maxBCG cluster sample. Solid 
line represents the best fitting model in case of three free pa- 
rameters: feefli 9i a-nd assumes "jackknife" error model. Long 
dashed line shows the same model without survey window con- 
volution applied and short dashed line corresponds to the linear 
model spectrum. Lower panel: the same as above, with the broad- 
band shape of the spectrum removed by dividing with a "smooth" 
model spectrum without BAO. Here we do not show linear model 
as it falls way out of this small plotted region. 



of maxBCG), n(r) is the radial selection function, and n is 
the area of the sky in steradians covered by the survey. 

In addition to the FKP errors we have also estimated 
the power spectrum variance via the "jackknife" method (see 
e.g. Lupton 1993). For this purpose we have divided the sur- 
vey volume into 5x5 angular regions (see Fig. \^ along with 
3 redshift bins, making the total of 75 regions, each contain- 
ing on average 168 clusters. Now the power spectrum is cal- 
culated 75 times with each time one of the regions omitted. 
The covariance of the power spectrum can now be estimated 
via the scatter of the power spectrum measurements: 

cov(P,Pj) = ^.^(p(")-a)(/^<"'-P,),(4) 

n = l 

1 ^ 

n=l 

Here P^"^ represents the power spectrum measurement for 
the i-th power spectrum bin in case the n-th survey subvol- 
ume was excluded from the calculations. All the 75 spectra 
along with the inferred "jackknife" error-bars for each of the 
spectrum bins i, APi — cov(Pi, Pi), are shown in Fig. |6l 
"Jackknife" and FKP error-bars are compared in Fig. [7] 
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Figure 5. Same as Fig. |4l with the only difference that here we 
have two free parameters: fe(,ff i 9i a^nd the value for u is fixed to 

ao/i-^Mpc. 
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Figure 6. Piecewise linear lines show 75 "jackknife" spectra and 
filled circles with error-bars represent the corresponding errors 
derived from the variability of these spectra. 

where for clarity we have removed the smooth broad-band 
component of the spectrum. We see that both errors are in 
good agreement. On the largest scales the "jackknife" errors 
probably slightly underestimate the real errors, whereas on 
smaller scales they might be more appropriate as they do not 
rely on the Gaussianity assumption. However, this small dif- 
ference does not influence our subsequent analysis, i.e. both 
error estimates lead to fully consistent results. 
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Figure 7. Comparison of the FKP and "jackknife" error esti- 
mates. 



Before we can fit model spectra to the data there are 
several effects one has to take into account: 

• Photometric redshift errors lead to significant damping 
of the spectrum; 

• Our observed spectrum is a pseudospectrum and thus 
before fitting one has to convolve model spectra with the 
survey window; 

• On smaller scales nonlinear effects become noticeable 
and one needs a model to account for these; 

• And finally, even without photometric redshift errors 
there are spectral modifications caused by the redshift-space 
distortions. 

Our model for the relation between the linear spectrum 
Piin(A:) and the observed one Pobs(fc) that tries to accom- 
modate all the four effects listed above reads as: 



s(fc) = J dk' k'^P{k')K{k',k) . 



where 



P(fc) = &cff(i + gfc^)/(fc)Piin(fc). 



(6) 



(7) 



Here b^s is the cluster bias parameter that also incorpo- 
rates the boost of the power at large scales due to the linear 
redshift-space distortions. Smaller scale redshift distortions 
are fully eliminated if the cluster detection algorithm is able 
to perform a perfect finger-of-God compression. However, in 
reality the compression is not ideal and thus there should be 
some extra uncertainties for the central redshift of the galaxy 
cluster. This can be modelled as an additional smoothing ef- 
fect in addition to the photometric redshift errors of the indi- 
vidual galaxies. The factor {1 + qk^^^) is our simplistic model 
to treat effects due to nonlinear evolution. This fo rm is very 
similar to the one suggested in ICole et all (|2005l ). however 
m iHutsil (|2006d ') we found that thepower law index 3/2 pro- 
vides better fit to the Halo Mode lj results compared to the 
value 2 used in lCole et all l|2005l ). The function f{k) in Eq. 



^ For a comprehensive review on Halo Model see lCoorav &: ShethI 
||2002|') . 



(O models the damping of the spectrum due to photometric 
redshift (photo-z) errors. Under the flat sky approximation 
it is easy to derive the analytic form for f{k) which reads: 



2ak 



erf((Tk) . 



(8) 



dt. This 



where the error function erf(x) = 
result assumes that photo-z errors follow Gaussian distribu- 
tion with dispersion Sz and a is the corresponding spatial 
smoothing scale, i.e. a = -^Sz . And finally, Eq. models 
the convolving effect of the survey window. We find that the 
mode coupling kernels K(k' , k) can be well described by the 
following analytic fit: 

K{k', k) = K{k, k') = ^ [g{k + k') - g{k - k')] , (9) 



where 



arctan 



kk' 



&4 



(10) 



and the normalization constant C is derived by demanding 
that / K{k, k')k''^ dk' = 1, giving: 

1 



C = 



(11) 



In the case of the maxBCG survey geometry the best fitting 
parameters a and b are found to be: a = 0.0044 an d b = 
0.0045. For the motivation of this analytic form see iHiitsil 
(|2006al ). 

In this paper we keep the background cosmology fixed 
to the "concordance" fiat ACDM model with Q,m ~ 0.27, 
Slfc = 0.045, and h = 0.7. The bias parameters quoted below 
are measured with respect to the model with erg = 0.85. Our 
decision not to fit for the cosmological parameters is based 
on two main reasons: 



(i) Photo-z uncertainties that are given in lKoester et al.l 
(|2007h are tested only for the case of the brightest cluster 
members that have spectroscopic redshifts available. The 
errors for the whole galaxy population in clusters might 
differ, and be possi bly larger than the numbers given in 
iKoester et al] l|2007l ): 

(ii) Uncertainties in the cluster finding algorithm which 
are hard to quantify and which might lead to the additional 
smoothing/overmerging, and thus to the extra uncertainties 
in derived cluster redshifts. 

Any untreated systematics in the determination of the ra- 
dial smoothing scale a will start to interfere with the cosmo- 
logical parameter estimation, resulting in biased estimates. 
Due to the uncertainties listed above we first treat the ra- 
dial smoothing scale cr as a free parameter. Thus we are left 
with three free parameters: 6cff , a, q. We also investigate 
separately the case where the radial smoothin g scale is fixed 
at c/H q ■ 5z ~ 30ft~^Mpc as suggested in K oester et al.l 
(|2007l ). Some of the results of these calculations are pre- 
sented in Figs. |4] and O and in Tables [2] and [S] To calculate 
theoretical model spectra we have used fitting formulae fo r 
the transfer functions presented in lEisenstein fc Hul (|l998l ). 
In the lower panels of Figs. |4] and [5] we have divided the 
spectra by the "smooth " model spectrum with acoustic os- 
cillations removed (see lEisenstein fc Hul [l99^ for details). 
For fitting purposes we have used Leven berg-Marqua r dt y^ - 
minimization technique as described in IPress et all (|1992| ) 
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Table 2. Best fitting model parameters and values along with the inferred confidence levels for BAO detection for the case with three 
free parameters: feeffi <?i 



P{k) 


3 free parameters: fe^g, 


g, (t; 16 — 


3 = 13 dof; 


expected — 13.0 ± 5.1 


error model 




(J 


•^wiggly 


Y'^ BAO detection 

^smooth 


FKP 


11.3 ±5.2 13.9 ±2.1 


64 ±30 


7.48 


10.9 1.8cr 


"Jackknife" 


8.7 ±2.5 14.1 ±2.7 


50± 16 


9.51 


14.2 2.2cr 



Table 3. Analog of Table [2] for the case with two free parameters: bcft i <?i a^nd the value for a is fixed to 30/i ^ Mpc. 



P(fc) 


u = 30; 2 free parameters: 




- 2 = 14 dof; 


expected — 14.0 ± 5.3 


error model 




-''■wiggly 


^smooth 


BAO detection 


FKP 


5.70 ±0.32 11.7 ±1.9 


11.8 


14.5 


1.60- 


"Jackknife" 


5.76 ±0.37 11.9 ±2.3 


12.5 


15.6 


1.8(T 



with the additional modification to allow for the nondiago- 
nal data covariance matrices. We perform the fitting for both 
error models: (i) FKP errors with diagonal covariance ma- 
trix, (ii) full "jackknife" covariance matrix. The solid lines 
in Figs. 2] and [S] represent the best fitting "wiggly" mod- 
els for the case of "jackknife" errors. In these figures we 
also demonstrate the damping effect of the window convolu- 
tion. The corresponding linear spectra are shown with short 
dashed lines. One can see that the shape of the measured 
spectrum deviates strongly from the linear model spectrum, 
which is mainly driven by the photo-z errors. The best fitting 
parameter values along with the corresponding values are 
given in Tables [2j [3] There we also present the value for 
the best fitting "smooth" model. In general we can see that 
the results for the x^-statistic agree well with expectations 
13.0±5.1 and 14.0±5.3 for 13 and 14 degrees of freedom, rele- 
vant for Table[2]and|31 respectively. Comparing the x^ values 
for the best fitting "smooth" and "wiggly" models one can 
assess the confidence level for the BAO detection, the results 
of which are presented in the last columns of Table [5] and 
[31 We note that using the full "jackknife" covariance matrix 
we obtain higher confidence for the detection of the oscilla- 
tory features in the spectrum. This is easy to understand; 
our first error model, where we assume diagonal covariance 
matrix and FKP errors is certainly not realistic. Survey win- 
dow, and more importantly nonlinear mode-mode coupling 
induces significant correlations between neighboring power 
spectrum bins. Stronger correlations amongst the bins make 
the spectrum "harder to distort" , and thus any surviving os- 
cillatory features gain higher statistical weight. It is interest- 
ing to note that in case we allow ct to be a free parameter the 
best fitting values are certainly larger than 30 Mpc. Also 
the bias parameter 6eff has a noticeable increase compared 
to the (T = 30 Mpc case. In fact there is a strong degen- 
eracy between &cff and a: if one chooses to increase 6cff , one 
can still obtain a good fit to the data by also increasing a by 
a suitable amount (see Fig. [8]). The fact that the case with 
three free parameters favors higher value for g (and thus also 
for 6eff ) is mainly driven by the relatively large amplitudes of 
the two first power spectrum bins. On the other hand, if one 
assumes that a = 30 Mpc is a reasonably correct value 
then one is left with some discrepancy on the largest scales 
between the model and data (see Fig. [5]). It is interesting 
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Figure 8. Error ellipses for the case with three free parameters. 
The solid and dashed lines correspond to the "jackknife" and 
FKP error models, respectively. The crosses mark the best fitting 
parameter values. Note the strong degeneracy between fe^g and 



that some large-scale "extra p ower" has also been detecte d 
by several other a uthors e.g. iPadmanabhan et al.l (|2006D ; 
iBlake et all \2m% who analyzed SDSS photometric LRG 
sampleQ Of course these mild discrepancies on the largest 
scales might be due to inaccurate treatment of the survey 
selection effects. Taking into account several uncertainties 
concerning the appropriate value for ct, and also uncertain- 
ties concerning the survey selection effects, we are here un- 
fortunately unable to draw any firm conclusions about the 
existence of the excess large-scale power. This issue certainly 
deserves a dedicated study and currently the best available 
sample for this purpose is probably the SDSS LRG spectro- 
scopic sample. 

However, one extra exercise we can perform here is to 
try to investigate how well the obtained values for the bias 
parameter 6eff agree with the model expectations. If one as- 



* See also Fig. 4 in lHutsilUoOGal where SDSS spectroscopic LRG 
sample reveals some mild excess of power on the largest scales. 
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0.1 0.15 
k [hMpc-^] 

Figure 9. Stability test of the power spectrum measurement with 
respect to changes in the radial selection function. The solid line 
with error-bars in the inset of the figure shows the measured 
power spectrum in case the radial selection corresponds to the 
volume-limited survey. The other solid line displays our "default" 
spectrum measurement. To increase the visibility of differences 
we have divided by the "default" spectrum in the main panel of 
the figure. 



sumes a flat ACDM model with the parameters as given 
above (i.e. i},n = 0.27, as — 0.85) then for the mass-limited 
cluster survey covering the same redshift range and sky area 
as maxBCG sample, and giving the same number of clusters, 
one infers a lower mass limit of Miow = 7.1 x 10^^ H'^Mq. 
We h ave calculated Miow usin g the Sheth-Tormen mass func- 
tion (jSheth fc Tormcn'' l999f ) averaged over the past light- 
cone (see , Hiitsi, 2006b for details). The corresponding light- 
cone averaged effective bias parameter including increase 
of po wer due to la rge-scale redshift-space distortions (see 
again iHiitsil I2OO6IJ 1 with respect to the z = linear spec- 
trum turns out to be foeff — 2.8 (b^g ~ 7.8)0 Comparing 
this value of b^g ~ 8 with the ones given in Tables [2] and 
[3] one sees that a = 30/i~^Mpc case has somewhat lower 
bias value, whereas the case with free a tends to have val- 
ues for fo^g, which seem to agree reasonably well (especially 
true for the case with "jackknife" errors). Thus this simple 
analysis would favor the "extra smoothing" hypothesis over 
the "extra power" one. 

As a final part of our power spectrum analysis we test 
the stability of our results against possible uncertainties in 



^ In reality of course for the mass-limited sample the number 
density drops with redshift, while maxBCG catalog is close to 
being volume-limited. We have calculated that to achieve a co- 
moving number density of 3.2 X 10~^ Mpc""^, as relevant for the 
maxBCG catalog, the lower mass limits for redshifts 0.1 and 0.3 
should be 7.8 X 10^^ h'^MQ and 6.6 X IQl^ h'^MQ, respectively. 
The effective bias parameters with respect to the z = linear 
matter power spectrum for these redshifts and lower mass limits 
are 6Jff°-^ ~ 2.75 and fe^f ~ 2.78, i.e. basically independent of 
redshift, and also agreeing very well with the light-cone averaged 
value of 2.8 quoted in the main text. This approximate constancy 
of the clustering amplitude at different redshifts for the objects 
with fixed comoving number density is a well known result. 




0.1 0.15 
k [hMpc-^] 

Figure 10. Analog of Fig.[9]for the stability test with respect to 
the variations in the angular mask. 



the survey selection function. For this purpose we have re- 
peated power spectrum calculations with modified radial 
and angular selections. In the inset of Fig. |5] we show the 
change induced in the power spectrum measurement if one 
replaces the smooth cubic spline fitted radial selection with 
the one corresponding to the volume-limited sample (see also 
Fig. [3} . For clarity the main panel of this figure shows the 
spectrum divided with our "default" spectrum. We see that 
at almost all the scales, except for the very largest ones, 
both spectra agree very well with each other. However, this 
exercise shows that the largest scales are very vulnerable 
with respect to the unaccounted slow trends in the survey 
selection. The results of the similar exercise for the angular 
selection function are shown in Fig. 1101 Here the agreement 
between the power spectrum measurements using the "min- 
imal" and "maximal" angular masks (see Section [2] for ex- 
planation) is almost ideal. Thus we can conclude that our 
results seem to be robust with respect to the uncertainties in 
the survey selection, which is especially true for scales with 
wavenumbers k > 0.03 ^Mpc~^. 



4 CONCLUSIONS AND DISCUSSION 

We have calculated the red shift-space pow er sp ectrum of the 
maxBCG cluster sample (Ko_ester_et_al. 20071 ) - currently 
the largest cluster catalog in existence. After correcting for 
the radial smoothing caused by the photometric redshift er- 
rors, treating the convolution with a survey window, and in- 
troducing a simple model for the nonlinear effects we show 
that "concordance" ACDM model is capable of providing 
very good fit to the data. Moreover, we are able to find 
weak evidence for the acoustic features in the spatial clus- 
tering pattern of the maxBCG sample. If radial smoothing 
scale (J, effective bias parameter b^s, and nonlinear distor- 
tion parameter q are treated as completely free parameters 
the model with acoustic oscillations is favored by 2.2cr over 
the corresponding "smooth" model without any oscillatory 
behavior. The evidence for BAO weakens somewhat (I.80"), 
and we are also left with some "extra power" on large scales. 
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if we fix the smoothing scale cr to 30 /i ^ Mpc as suggested 
by the photometric redshift errors 5z = 0.01 estimated in 
iKoester et al] (|200'it . We note that these photometric er- 
rors actually apply only for the very brightest cluster mem- 
bers and in reality the scatter for the whole sample might 
be larger. Also the cluster finding algorithm itself might in- 
troduce additional smoothing/overmerging along the line of 
sight. Due to these uncertainties we have decided not to 
carry out any cosmological parameter study in this paper, as 
any untreated systematics will be immediately propagated 
to the parameter estimates, leading to biased results. We 
have estimated power spectrum errors using two different 
methods: (i) FKP errors, (ii) "jackknife" errors, obtained 
by dividing the survey into 5x5x3 chunks containing 
equal number of galaxies. Both of these errors turn out to 
agree very well with each other. We have also shown the sta- 
bility of our power spectrum measurement with respect to 
the possible uncertainties in the radial and angular selection 
function. 

The detectability of BAO with ~ 10'' galaxy clusters 
is very encouraging result keeping in mind the future clus- 
ter surveys, such as the ones based o n the measurement of 
the t hermal Sunyaev-Zeldovich effect jSunvaev fc ZeldovichI 
Il972h . e.g. Spt0, PLANCiMor the proposed 100,000- 
cluster X-ray survey eROSIT43- In reality the main target 
of these cluster surveys is to map the evolution of the cluster 
number density as a function of redshift, since this quantity 
is very sensitive to the amplitude and to the growth rate of 
the density perturbations, and as such it provides a powerful 
tool to constrain the properties of DE. The possibility to also 
find the clustering power spectrum can be seen as a useful 
byproduct of these surveys that comes essentially "for free" . 
However, both cluster number count study and spatial clus- 
tering analysis needs as an input estimates for the cluster 
redshifts, with the former probably doing relatively well with 
rather poor redshift estimates (e.g. photometric redshifts). 
This is a rather complicated task since one arguably needs of 
order 10 galaxy redshifts per cluster to reliably infer cluster 
redshift. In that respect for the measurement of BAO it is 
certainly less costly to sample the cosmic density field using 
e.g. LRGs or blue emission line galaxies. The latter type of 
objects are targets for th e currently ongoing WiggleZ project 
ijClazebrook et aLll2007^ at AAO that attempts to measure 
redshifts for 400, 000 objects over ~ 1000 deg^ and cover 
the redshift range 0.5 < z < 1. Ther e is also a Wfmos in- 
strument (see e.g. iBassett et ahlfioosl ) construction planned 
for the Gemini and Subaru observatories that will hopefully 
be completed in 2012. This new multi-object spectrograph 
will be able to measure the spectra of ~ 5000 objects at the 
same time. There are plans to perform a wide field (~ 2000 
deg^) redshift survey giving spectra for ^ 2 x 10*' galax- 
ies up to redshifts of z ~ 1.3 together with a narrower 
(~ 200 deg'^) and deeper (z ~ 2 — 3) survey with a yield 
of ~ 5 x 10^ galaxies. For the more distant future (around 
2020) one would expect a superb measurement of the mat- 
ter power spectrum using ~ 10® galaxy redshifts obtainable 
via the 21 cm measurements by the Square Kilometre Array 



(SKAjB (e.g. iBlake et al.ll2004 lAbdaUa fc Rawhndlioosl ). 

However, before Wfmos and Ska become available several 
wide area imaging surveys, such as Pan- STARR^3 or Dark 
Energy SurveM^j (Des), will be performed. Although these 
photometric surveys lack accurate redshift information, the 
huge number of objects detectable over large sky areas signif- 
icantly compensate this shortcoming, allowi ng us to obtain 
a highly competitive measurement of BAO ([Blake fc Bridle! 
I2005h . 
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